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I [OIII]5007A using data from SDSS DR7. A special attention is paid to a subclass of 
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ly-^ . narrow broad permitted emission lines with FWHM < 2000 km/s. A new approach 
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_Cy _' BLSyl, but the average values are within errors. 
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1. INTRODUCTION 

It is not yet completely clear how do active nuclei (AGN) and their host galaxies evolve. 
During the recent years there is a growing tendency to pay a close attention to secular 
processes. Such processes are quite slow and are caused by different types of internal insta- 
bilities in contrast to mergings. For example, [1, 2] showed that there is a high fraction of 
systems with disk-like morphology of host galaxies among high-redshift quasars and with a 
high amount of gas in them. 

According to [3], at redshifts z ~ 2, many of the galaxies with "clumpy" morphology that 
thought to be the results of mergings show regular velocity maps similar to the normal disk 
galaxies. The observed morphology is a result of a very active star formation caused by gas 
infall along the large-scale filaments. 

It is believed that the growth of a central SuperMassive Black Hole (SMBH) must be 
connected with the growth of the host galaxy or with some of its components — bulge, 
halo etc. Many researchers reported a very tight correlation between mass of the spherical 
components of galaxies and masses of their black holes (see, for example, the classical paper 
[4]) It is also accepted that there is a correlation between the mass of central black hole and 
the velocity dispersion of stars in bulge: Mbh ^ cr* , where n > 1 (see, for example, [5, 6]). 
A "Spherical component" is the whole galaxy for elliptical systems and the bulge for disk 
galaxies. 

Nevertheless, many galaxies have the so-called "pseudobulges". They have lower S/T ra- 
tios compared to classical bulges, low Sersic indexes {ni, < 2) and also high angular momenta. 
Such bulges are dynamically more similar to disks of galaxies than to spherical non-rotating 
systems. It is believed that pseudobulges are formed via secular evolution, which in turn can 
be triggered by different internal processes or by outer influences such as flyby of a satellite 
galaxy or a minor merging [7]. It is very important to stress that these bulges do not follow 
the well-known correlations Mbh - M^uige and Mbh - cr*- 

Narrow-Line Seyfert galaxies are good representatives of such type of systems. They were 
first recognized as a distinct type of active galaxies in [8] as Seyfert galaxies having unusually 
narrow broad permitted hues in their spectra (FWHM < 2000 km/s). Ordinary Seyfert 
galaxies in context of this classification are called Broad-Line Seyfert galaxies (BLS). NLS 
show many other peculiar properties. Host galaxies of NLS are, on average, of later types 
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than those of BLS. The average Hubble Types of host galaxies are (T) — 3.0 and (T) — 1.0 
for NLS and BLS, respectively [9]. Numerical encoding of the Hubble types is according to 
the Third Reference Catalogue (RC3) [10]. T — 1 corresponds to Sa galaxies, T — 2 — Sab, 
T = 3 — Sb and so forth. 

Their host galaxies always have pseudobulges judging on NLS galaxies with available 
photometry. On the other hand, the host galaxies of BLS may have classical bulges as well 
[11]. 

It is now widely accepted that NLS accrete at high Eddington ratios, close to unity and 
that their black holes are undermassive when compared with those of BLS [12, 13]. Properties 
of NLS in 7", x-ray and radio deserve a special attention. Many researchers have dedicated 
their studies to this problem. The main result of these studies is the recognition of the 
presence of relativistic jets in NLS. The very existence of jets in late- type spiral galaxies with 
low mass black holes is an intriguing discovery. Some authors supposed that high angular 
momentum of supermassive black holes is a factor that leads to the launch of relativistic jets 
in such systems. Such a rapid rotation is acquired by a black hole via accretion of gas with 
a high angular momentum, transported to the galactic center by secular processes. 

2. DATA PROCESSING 

For our study we use data from the SDSS DR 7 (Sloan Digital Sky Survey Data Release 
7) [14]. It needs to be mentioned that in our study we do not use diagnostic diagrams [15, 16] 
because there is no necessity in it. This work is dedicated to AGN type I with widths of the 
broad permitted lines FWHM > 1200 km/s. Special methods of distinguishing there objects 
from SF (Star Forming) and TO (Transition Objects) are not required. 

We use Ho; line instead of H/3 to classify AGN into NLS and BLS. It is known that H/3 line 
has a prominent narrow component. This is why it is impossible to use width estimations of 
this hue based on a single Gaussian fit from SDSS catalogue to classify AGN. In figure la the 
FWHM(H^)sDss is plotted versus FWRM(Rl3)broad- FWHM(H^)sdss is the H^ width from 
SDSS survey, FWliM{E.p)broad is the FWHM of the broad component of H^ line from [17], 
where authors thoroughly fitted the spectra. As can be seen from this figure, it is impossible 
to use widths of H/3 line from SDSS. 

Let us now direct our attention to the Ha line. It does not show such a prominent narrow 
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component. In figure lb FWHM(Hq;)sdss is plotted versus FWHM(H^)broad- It is necessary 
to exclude objects which have at least one line narrower than the threshold value 1200 km/s. 
Objects without Broad Line Region are AGN type 2 and are not subject of this study. As 

can be seen in fig. lb there is a tight correlation FWHM(Hq;)sdss = b ■ FWB.M{E./3) broad + o, 
with the following parameters: 

a = 0.499 ±0.011, 6 = 906 ± 46 . 

The deviation of a coefficient from unity can be easily explained. The excitation energies 
of Ha and H/i lines are different. That is why Ha line will be emitted at larger distances 
from the active nucleus and will have lesser FWHM. For this study the crucial moment 
is the presence of a correlation which leads to a possibihty to classify AGN into BL AGN 
(Broad-Line AGN) and NL AGN (Narrow-Line AGN) using the data on Ha widths. 

The task of calculation of the intrinsic AGN luminosity is not among the ultimately solved 
ones. All approaches described in literature can be divided in two groups: based on emission 
lines and continuum. AGN luminosity can be estimated using continuum observations only 
when host galaxy contribution is negligible. For luminous quasars it is the case. For exam- 
ple, such an approach was used in [18]. In order to measure luminosities of fainter AGN is 
it necessary to use photometric fitting. In different papers were used different sets of com- 
ponents: Point Spread Function (PSF) + Sersic profile [2], PSF + two Sersic profiles with 
different rib [19], PSF + de-Vaucouleur profile + exponential disk [20]. But this approach 
has several significant drawbacks. At first, many AGN lack photometric observations of suf- 
ficient quality to separate photometric components. This is a significant problem for SDSS 
survey because resolution of it is not very high. Secondly, proper accounting for orientation 
effects and internal absorption is a difficult and often an impossible task. 

The use of emission fines to estimate AGN luminosity has some significant advantages. 
In this study we utilize forbidden oxygen line [OIII]A5007A. In some papers it is stated that 
the contribution from star formation to this line is negligible (for example, [21]). This line is 
formed in the Narrow-Line Region which is more distant from the nucleus than BLS and has 
a double-cone geometry. This implies that observed luminosity in [OIII] does not depend on 
orientation [21]. 

In this study we do not perform the spectral fitting. It is necessary to check if it is 
possible to use line intensity data from SDSS. In order to do so we plot lg/[oiii],SDSS versus 
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lg-?"[oiii],[i7] (see fig. 2). lg/[oiii],SDSS is an [OIII] intensity from SDSS survey, Ig /[oiii], [17] 
is an [OIII] intensity obtained in [17] by accurate spectral fits. The relation is as follows: 
Ig /[OIII], SDSS = lg-^[oiii], [17] + a- , where a = —0.0053 ± 0.0027. Due to the strength of 
the correlation [OIII] intensities obtained by the automatic SDSS pipeline can be used to 
calculate AGN luminosities. 

We select total 9020 BL AGN objects from SDSS DR7 with detectable [OIII] and Ha 
emission having FWHM(Hq;) > 1200 km/s. We divide this sample into two subsamples, 
NLS and BLS. There are 2082 and 6938 objects in these subsamples, respectively. We set 
upper limit on the redshift equal to 0.18 because on higher redshifts it is impossible to apply 
our normalization algorithm (see details below). The lower limit on the distance was set to 
55 Mpc to avoid the influence of the Local Supercluster. In the redshift interval considered 
the spectral lines used in this study are within the interval of spectral sensitivity of the SDSS 
survey. 

In this study we used the following cosmological parameters in accordance with SDSS 
DR7: Qm = 0.279, = 0.721, h = 0.701. Results of other studies were converted to this 
cosmology. 



3. THE LUMINOSITY FUNCTION 

The are two main methods of obtaining the luminosity function. These are the maximum 

V 

likelihood method and the — — method. In the recent years due to substantial growth of 

''max 

the volumes of the catalogues there is a growing tendency to use a more simple, but having 

V 

significant advantages — — method. The main advantage is that this method is free from a 

*max 

V 

priori presumptions about the shape of the luminosity function. The — — method implies 

'max 

calculation for each object the maximum volume where it can be detected by the observer. 

The task of calculating the Vmax from emission fine observations is quite complex. Because 

there are several reasons why an object cannot be observed: visual magnitude is of limits of 

completeness of the survey, the decrease of S/N with distance, widening of the lines at larger 

V 

redshifts, pixel overfill in spectra of close luminous objects etc. In order to apply the — — 

<'max 

method one must create model spectra for each object, simulate the noise and then calculate 
maximum volume from which the object under consideration can be observed. Another 
drawback is that it does not take into account variations of the Large Scale structure. 
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In order to use our method it is necessary to obtain a function of average probability of 
object detection p{dc), where dc is the comoving density. Hence, we can derive average value 
of Vmax for objccts in the considered luminosity bin 



r^c, max 
" dr.. m.in. 



where dunlin and dc,max are the lower and upper distance limits, respectively. The value of 
the luminosity function in the considered bin is: 



'norm 



(Vmax) d Ig -L[OIII] ' 

where Nnorm is the amount of galaxies after normalization. The (p{L) is expressed in the 
units of Mpc~^ (lgL[oiii])~^ i.e. the amount of objects per unit of volume per unit of the 
logarithm of the luminosity. 

The total error of estimation consists of the Poisson error (TN„^^,poiss — V^obs and the 
error in estimating (vmax) {'^(vmax)) and is equal to: 



^max ) ^ 



cilgL[om] V V {'"max) J \ {v 
where ^ is a factor corresponding to the normalization, so N„orm — d^obs- In this study 
all the luminosities are expressed in units of L/Lq, where Lq — 3.844 x 10^^ erg/s. As is 
said above, the existing methods do not take into account the variations of the Large-Scale 
structure of the Universe. In order to account for this effect we use the following approach. 

The luminosity function of inactive galaxies for the Local Universe is well known and was 
obtained in many independent studies. Let us define {pgai) as the average density of galaxies 
calculated from the luminosity function and Pgai as an observed density of galaxies in the 
considered volume. The amount of AGN satisfying some predefined conditions is Nagn- 
This allows us to obtain the amount of AGN in the considered volume normalized to the 
average density of galaxies in the Local Universe: 

Nagn, norm — Nagn . 

Pgai 

In order to proceed the normalization let us create a sample of galaxies limited in magnitude 
by rrimin — 14.5 and rrimax — 17.6. Within these hmits the SDSS survey is beheved to 
be complete. Let us define Lghg as observed total luminosity of inactive galaxies in the 
considered volume. If Lj is a luminosity of i-th galaxy, then L^bs = S^li L^. 
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The luminosity function is usually fitted by 

4){L)dL = (j)*{L/L*)"e^p{-L/L*)d{L/L*) 

or in the absolute magnitudes [22] 

(j){M)dM = (O.41nlO)0*lO°-^(*^*-^)(^+")exp (-10°-^(^*-^))(iM . 

In [23] the following parameters of the local luminosity function in r band of the SDSS are 
obtained: 0* = 0.0090 ± 0.0007, M* -5\gh = -20.73 ± 0.04, a = -1.23 ± 0.02 Using these 
parameters it is possible to obtain the calculated total luminosity of galaxies per unit of 
volume: 

Lsch= [ ' L'(j>{L')dL. 

J Lx 

Integration in this formula should have upper and lower limits L^ax and L^m oi' ^max and 
Mmin- It is because for faint galaxies the luminosity function is not well defined and when 
high value of L^^ax is chosen the result of integration drastically depends on the selected 
parameters of the Schechter function. In this study we used M^ax — —17, = —23 in 

the r band. We also use these hmits when we calculate the Lohs- 

It is now possible to obtain the relation of the average density of galaxies in Local Universe 
to observed density of galaxies in considered volume: 

{pgal) LschVp 
= — = K, 

Pgal ^ohs 

where V is the volume of the Universe within distance limits dc,min,dc,max S'Hd // is the 
coverage of the survey. 

Of course, this relation, is correct only when parameters of the Schechter function are 
constant. For galaxies in the Local Universe it is the case. We consider the normalization 
possible if Lsch/ Lsch,tot is larger than threshold value 0.1, where Lsch,tot is the total luminosity 
of all galaxies, derived from Schechter function. It quite interesting that when this approach 
is used it is not required to know the coverage area of the survey. Since when concentrations 
are used we have: 

LschVfi 1 Lsch 

n AG N, norm — AGN — ^ TT — ^^AGNJ ■ 

J^obs P ^ J^obs 

In means that Lobs already contains the survey coverage in it. Let us now check the reliability 
of the described algorithm. In order to do so we apply it to a sample of galaxies with known 
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absolute magnitudes. As we have already mentioned above, the SDSS survey is complete in 
the magnitude interval 14.5 < < 17.6 [14]. The dependence of galaxy concentration on 
z in Ml < M < M2 interval can be expressed in the following form: p{z) = p^fgs{z)fc{z), 
where fgs{z) is a term that reflects variations of the Large Scale structure and /c is the 
completeness of the survey in selected interval on certain z. After the normalization we 
obtain: p{z) — p*fc{z)- If the normalization algorithm is working properly, then a plateau 
followed by a rapid decline when visual magnitude reaches the completeness limit (see fig. 3a) 
of SDSS survey should be seen. This is the exact description of what can be seen in fig. 3b. 

The luminosity function of inactive galaxies can be obtained directly. The concentration 
of galaxies in absolute magnitude interval Mi < M < M2 is equal to the plateau of the 
p{z) function. The situation with AGN is far more complex. The are many reasons why an 
object will not be in the AGN sample. These reasons can be divided in two groups. The first 
group consists of factors that make observation of close objects impossible (too bright visual 
magnitudes, very bright emission lines that leads to overfill of pixels in their spectra etc.). 
The probability function Pinc{dc) will increase with distance. Let us search it in the form 
exp {—bi/dl), where bi is a free parameter. The second group consists of factors that cause 
the decrease of detection probability with distance due to decrease of the observed flux. The 
corresponding function p^ec will decrease with distance to the object. Let us search it in the 
from of exp {—b2dl), where 62 is a free parameter. 

Now we will briefly address the issue of misclassification. The characteristic that defines 
the NLS and BLS populations is the width of broad permitted lines, in our case Ha. As was 
noted in [24] when the signal to noise ratio is low FWHM of emission lines is systematically 
underestimated. That is why there will be a tendency to misclassify a genuine BLS due to 
the underestimation of FWHM with increasing redshift. Some of the genuine NLS will be 
misclassified because their measured FWHM will actually be below the threshold between 
Syl and Sy2. Also some of objects which are actually BLS will be classified as NLS at higher 
redshifts. All these effects will contribute to considered probability functions. 

For each luminosity bin L[oui] ^ [Li, i^i+i] the normahzed AGN density must be fitted 
with the function pagn = aPdecPinc- The corresponding average value of v^ax is: 



Because the observed values of pagn{z) are actually small, especially for luminous objects 



max 
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we use cumulative function: 

fobsidc) — NAGN{dc,i < dc) , 

where NAGN{dc,i < dc) is an amount of AGN in the considered luminosity bin within distance 
limits dc,min < dc,i < dc- For fitting we use the following function: 

ffit{dc,i)^ 47radle(~i)e^-''''"^d{dc). 

It should be noted that the main result of the approximation is the shape of the function 
PAGN{dc)- The value (T{vmax) derived from the error of the approximation, i.e. from (7(0), 
(t(6i) and (7(62) (see below). 

We use neither a fixed step in luminosity nor in the distance. The bin widths were 
selected manually because the luminosity function varies drastically with iv[oni]- Bin widths 
and amounts of objects in every bin are listed in table 1. A fixed step in distance makes 
it impossible to use our normalization algorithm at low distances. The volume in each bin 
depends on the distance and the amount of objects in close bins is actually low. Instead of 
it we use a fixed step in volume 5Vc. For the i-th bin we have the distance limits: 

Finally, we obtain luminosity functions for three samples: all Syl, NLSyl, BLSyl. 

All calculations are done using programs written in PYTHON language by the author of 
this paper. 

4. A COMPARISON OF THE OBTAINED LUMINOSITY FUNCTION WITH ONES 

FROM THE LITERATURE 

The luminosity functions of NLSyl, BLSyl and all Syl are shown in fig. 4. The results 
obtained are also listed in table 1. The luminosity function from [21] and [25] are showed 
for comparison. In [21] LFs were obtained for Syl and Sy2, in [25] only for Sy2. As can 
be seen in fig. 4, the luminosity functions from [21, 25] are in a good agreement with each 
other. Our luminosity function is significantly lower. 

But, the luminosity functions from many different studies show disagreements with [21]. 
In order to compare our results with as many as possible LFs calculated by different authors 



10 



we convert our LF from [OIII]A5007A to the Johnson band B according to [25]: 

Ig = -0.38Mb - 0.42 , 

In order to convert luminosities from Ho; to Mb we use the relation from [26]: 

Mb = -2.1(lgLHa-42) - 20.1. 

It should be noted that these relations have significant dispersions. For example, it results 
in disagreement between luminosity functions converted to B band from Ha and [OIII] from 
the same study [21]. 

The result of conversion of luminosity functions is shown in fig. 5, where the following 
LFs are also plotted: 

1. Converted from [OIII] to Mb ([21, 25] and this study) 

2. Converted from [21, 26, 27] 

3. The quasar luminosity function at redshift 0.4 < 2; < 0.88 measured in B band [18]. 

There is a known erratum in [27], we used corrected data published by authors in [28]. 

As can be seen from these figs LFs obtained by different authors form two distinct groups 
(groups II and III in fig. 5). They consist of the results of [21, 25] and [26, 27], respectfully. 
LF of Syl obtained by us is between these contradicting results. It is important that our LF 
agrees well with the LF of quasars [18] on highest luminosities. 

It is remarkable that (f){NLS)/(j){BLS) is not constant. Instead, it is actually a function 
of the luminosity (fig. 4., bottom panel). It has a maximum between log {L[oiii]/Lq) — 6.8 
and 7.6. At larger luminosities it slightly decreases. The decrease towards lower luminosities 
is ~ 0.5dex to \og{Loiii/ Lq) — 5.5. 

A putative explanation can be proposed when the results of [29] are considered. The 
authors studied samples of NLS and BLS selected from SDSS DR3. The AGN luminosities 
were estimated using pgp and pgp. That is why only objects with L^gn ^ Lhg 
were included in the sample. The ones of the key results of their work are the distributions 
of log {L/Lgdd) S'Hd log (Mbh/Mq). As was stated in [29] and many other studies, the 
distributions of Mbh and L/Ledd differ significantly for NLS and BLS. 
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Due to the fact that Ig L = Ig (ledd = Ig (1-3 x 10=^«)+lg ^+lg y^, if Ig ^ 

L 

and Ig - — have normal distributions Ig L also has a normal distribution. 

Ledd 

The parameters of Ig L^; distributions are: 
for NLS: 

/X = 44.589 ± 0.052, 
a = 0.404 ± 0.043, 

for BLSyl: 

II = 44.535 ± 0.088, 
(7 = 0.560 ± 0.061, 

where // and a are the expected value and dispersion. The expected values are within the 
errors, but dispersions differ significantly. Thus, (l){NLS)/(f){BLS) ratio has a maximum 
at log (Lboi) — 44.64 which corresponds to log {L^oiii]/ Lq) — 7.4. The dispersion of NLS 
distribution is lower, therefore NLS occupy a narrower range in luminosity than BLS. 

We would hke to stress it again that the comparison with [29] is putative. In [29] Lagn 
was estimated using M^, sdss and sdsS) i-©- the distribution of Lagn in their sample is 
different from this study. 

We calculate a predicted soft x-ray(0.5-2kev) luminosity function based on the L^oiii] 
data. The luminosities in [OIII]A5007A are converted into the soft x-ray following [30, 31]. 
The result is plotted in fig. 6. The x-ray LF of Syl based on the ROSAT, XMM-Newton 
and CHANDRA data [32] is also shown in fig. 6. For a comparison the predictions based on 
the Ha luminosity functions from [26, 28] and [OIII] from [21] are also plotted in fig. 6. 
luminosities are converted to the bolometric ones following [27], then from the bolometric 
ones to Lo.5-2feei; following [30]. 

It should be stressed that our prediction corresponds to the actually observed x-ray lu- 
minosity function better than the LFs from other studies. 

5. CONCLUSIONS 

• We developed a method of evaluation of AGN luminosity function based on emission- 
line data accounting for variations of the density of galaxies due to the large-scale 
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structure. The [OIII]A5007A emission-line luminosity functions for NLS, BLS, Syl 
(NLS+BLS) are obtained. 

• In order to compare our results with luminosity functions from other studies we con- 
verted our LF to the B band. LF of Syl obtained by us is between contradicting 
results from different pubhcations. 

• At largest luminosities the LF of Syl agrees very well with the one of quasars [18]. 

• A prediction of soft x-ray luminosity function based on our emission-line LF is in an 
agreement with the observed x-ray LF. 

• (f){NLS) I (f){BLS) is not a constant, but a function of luminosity. It has a maximum be- 
tween log (L[o///]/Lq) = 6.8 and 7.6. This agrees with the results of [29]. The average 
luminosities of NLS and BLS are equal within errors, but widths of the distributions 
differ. NLS occupy a narrower range of the AGN luminosity. 
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Table 1. The luminosity functions of NLS, BLS and Syl(NLS + BLS). The units are as following. i/[oiii] 
is in units of Igi/L©, Bin (size of the luminosity bun) is in dex, luminosity function (j>{L) in 
Mpc~^(lgL[oiii])^) ^ is an amount of objects in every luminosity bin. 



NLSyl 



BLSyl 



Syl 



^[OIII] 



Bin 



4>{L) N 



^[OIII] 



Bin 



0(L) 



L 



[OIII] 
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6.875 
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7.625 
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8.25 
8.75 



0.5 
0.25 
0.25 
0.25 
0.25 
0.25 
0.25 
0.25 
0.25 
0.25 
0.5 
0.5 



-5.138 



+0. 



-5.150t[!; 
-5.1541°; 
-5.213tg; 
-5.399tg; 



-5.190^ 



-5.463 



-fO. 
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-5.51llg; 
-5.80;0; 
-6.24ll°; 
-6.76^ 
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10 
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'.14 
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5.4375 
5.5625 
5.6875 
5.8125 
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6.1875 
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6.4375 
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0.25 
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023 

020 
021 

1.021 
1.022 

021 
022 

023 
025 

026 
028 

028 
029 

024 
026 

025 
027 

023 
025 

030 
036 

023 
1.024 

1.049 
1.055 

029 
031 

064 
075 



299 
323 
388 
439 
469 
470 
443 
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307 
346 
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0.125 
0.125 
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Figure 1. The dependence of the FWHM of the H/3(a) and Ha(b) based on SDSS single Gaussian 
fitting on the FWHM of the H/3 from [17]. Grey dots indicate the objects without BLR for one of the 
lines, i.e. FWHM(H/3) < 1200 km/s or FWHM(Ha) < 1200 km/s (horizontal line). 




-16 -15 -14 -13 

log I[,7j (erg s"'cm"^) 



Figure 2. A comparison between integral intensity in [OIII] line obtained by the automatic pipeline 
of SDSS with the intensities taken from [17] obtained by a thorough spectral fitting. 
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Figure 3. (a) — The dependence of absolute magnitudes on the comoving distances for visible 
magnitudes corresponding to the completeness limits of SDSS. The horizontal dashed lines are the 
limiting absolute magnitudes used in this study for normalization. See details in the text, (b) — The 
dependence of normalized volume density on the comoving distance for different narrow intervals in 

the absolute magnitudes. 
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Figure 4. The luminosity functions of NLS, BLS and Syl(NLS + BLS) obtained in this study. For 
the sake of comparison LFs from [21] and [25] obtained in the same line [OIII[ A5007A are also plotted: 
1 - Syl ]21]. ^- Sy2 ]21]. 5- AGN type 2 ]25]. 4, 5, 6- NLS, BLS, Syl(NLS + BLS), respectively, 
obtained in this study. Group I consists of luminosity functions from the literature. Group II — 

luminosity functions obtained in this study. 
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Figure 5. The luminosity functions in the B band obtained by converting luminosities from [OIII] 
A5007A or Ha: 1, 2, 3 - NLSyl, BLSyl, NLSyl + BLSyl, respectively, obtained in this study. 4 
- AGN type 2 from [OIII] [25]. 5 - Syl from ]OIII] ]21]. 6 - Sy2 from ]OIII] ]21]. 7 - Syl + Sy2 
from Ha [21[. 8 - Syl from Ha [26[. 9 - Syl from Ha [28[. 10 - QSO, obtained in B band [18]. 
The luminosity functions are divided in four groups. Group I consists of LFs obtained in this study. 
Group II — LFs from [21, 25[. Group HI — LFs from [26, 27[. LFs in groups II and HI show a good 
agreement with each other. Group IV — the LF of quasars [18[. 
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Figure 6. The observed Syl x-ray Luminosity Function (0.5-2 kev) [32]. The predicted Syl x-ray 
LFs based on the luminosity functions from this study and the ones from the literature are also plotted; 
1 - x-ray 0.5-2 keV [32]. 2 - hY prediction based on [OIII] LF from this study. 3 - ]OIII] ]21]. ^ - 

Ha [26[. 5 - Ha [28[. 



